
s1 = simulate(m,d,1:T,'anticipate=',true,'nonlinear=',T,'maxiter=',1000);

s2 = simulate(m,d,1:T,'anticipate=',true,'nonlinear=',T,'maxiter=',1000, ...
    'sparseShocks=',true);
